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Abstract 


In this paper, we reconsider the single interval Renyi entropy of a free compact scalar 
on a torus. In this case, the contribution to the entropy could be decomposed into classical 
part and quantum part. The classical part includes the contribution from all the saddle 
points, while the quantum part is universal. After considering a different monodromy 
condition from the one in the literature, we re-evaluate the classical part of the Renyi 
entropy. Moreover, we expand the entropy in the low temperature limit and find the leading 
thermal correction term which is consistent with the universal behavior suggested in 15 


Furthermore we investigate the large interval behavior of the entanglement entropy and 
show that the universal relation between the entanglement entropy and thermal entropy 
holds in this case. 
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1 Introduction 


Entanglement entropy is an important notion in many-body quantum system Si- For a 
bipartite system, the entanglement entropy of the subsystem A is defined to be the Von 
Neumann entropy of the reduced density matrix 

See = -Tr pa log pa, (1-1) 

where the reduced density matrix pA = Tr^ p is obtained by smearing the degrees of freedom in 
B complementary to A. When the total system is in the vacuum p =| 0)(0 |, the entanglement 
entropy of A equals to the one of its complement 

See{A) = See{B). (1.2) 

At a finite temperature, however, such equality does not hold anymore. The entanglement 
entropy has been studied in various condensed matter systems [^, and also in the context of 
AdS/CFT correspondence (^[^. 

The computation of the entanglement entropy from its definition becomes a formidable 
task when the number of degrees of freedom in the system is huge. Especially in quantum field 
theory with infinite degrees of freedom, it is more convenient to compute the Renyi entropy 
via the Replica trick j^. The Renyi entropy is defined to be 

Sn = -^Trpd. (1.3) 

It is related to the entanglement entropy by the relation 


See = lim Sn 

n—>-1 


(1.4) 


if the analytic continuation of the limit is available. In field theory, the entropy is defined 
with respect to a spacial submanifold at a fixed time. In two dimensional spacetime, the 
submanifold could be interval(s). However, after Wick rotation, the Euclideanized field theory 
is defined in a complex plane and the Renyi entropy becomes 


Sn = - 


1 Zn 
n- 1^’ 


(1.5) 


in which is the partition function on a n-sheeted Riemann surface resulted from the pasting 
of n complex plane along the branch cuts (intervals). 

There has been a long history to study the entanglement entropy in quantum field theory. 
In the dimension higher than two, the entanglement entropy has been found to obey the area 
law (for a nice review, see [^). But we are only allowed to compute the entanglement entropy 
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analytically in very restricted situations, for example the one for sphere in free field theory. In 
(1 + 1) dimension, we can do better, especially for the field theory with conformal symmetry. 
For a 2D CFT on complex plane the Renyi entropy for one interval of length i is universal 
and only depends on the central charge 


n + 1 I 

= c— —log-. 
bn € 


( 1 . 6 ) 


For more complicated cases, for example the multi-interval at zero temperature or single 
interval on a circle at finite temperature, the entanglement entropy and Renyi entropy depend 
on the details of the CFT. In [9,10, the double interval Renyi entropies for free bosons and 


Ising model have been computed. In 11,12 the finite temperature Renyi entropy for free 


fermions has been discussed. In 13 , the finite temperature Renyi entropy for free bosons has 
been studied. Moreover, the treatment on the IT-functions in the partition functions of free 
boson has been improved in [14| . 

In this paper, we reconsider the Renyi entropy of free compact scalar on a torus. Our 
motivation is two-fold. First of all, in our study on the free noncompact scalar case [^, we 
noticed at least two novel features, originated from the continuous spectrum of the theory. One 
is that the leading thermal correction at low temperature in this case takes a form different 
from the one suggested in |15|. This is because that the noncompact scalar has a degenerate 


vacuum, while the universal thermal correction found in 15 was based on the assnmption 
that the CFT has a mass gap. The other one is that the universal relation between the large 
interval entanglement entropy and thermal entropy does not hold any more |14[|16[ . Since the 
noncompact free scalar could be taken as the large radius limit of the compact free scalar, it 
would be interesting to study the Renyi entropy of free compact scalar, which has a discrete 
spectrum. On the other hand, even though the Renyi entropy of free compact scalar has been 


computed in 13 , the detailed discussion on the low temperature or large interval expansion 
has not been worked out. Moreover, we find that the classical part of the partition function 
actually depends on a different monodromy condition from the one in [^ . With the corrected 
classical partition function, we compnte the Renyi entropy and do expansion in several limits, 
and rediscover the expected universal behaviors. 

In the next section, we revaluate the Renyi entropy for free compact boson. We consider 
a slightly different monodromy condition to read the classical part of the partition function. 
Then in Section 3, we discuss a low temperature limit, and expand the Renyi entropy with 
respect to We find the leading order is now consistent with the universal thermal 

correction suggested in [15] . In Section 4, we investigate the small interval and large interval 
limits of the Renyi entropy, and prove the universal relation between entanglement and thermal 
entropies. 
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2 Compact boson Renyi entropy 


For a free boson, the partition function on a Riemann surface can be decomposed into classical 
and quantum parts 

( 2 . 1 ) 


^ ^quantum^classical' 


The classical part gets contributions from all the saddle points carrying different monodromy 
conditions 

Zclassical = Y.^~^''^^^- ( 2 - 2 ) 

M 

For quantum correction, we need to consider the perturbations around the classical saddle 
point and evaluate their contribution to the partition function. In general, for different classical 
solution their quantum corrections are different, but in the case of free bosons the quantum 
correction is universal so the partition function can be decomposed into the classical and 


quantum part as in (2.1). 


In this section, we compute the compact scalar partition functions on n-sheeted torus 
connected by a branch cut. The free complex scalar is compactified on a square torus of radius 
R. It obeys the boundary condition 


X{e‘^'^^z,e = X{z, z) + 2TrR{mi + im 2 ), mi,m 2 ^Z. 


The quantum part of the partition function equals to [13, 14 

n—1 


z. 


1 


n, quantum 


n 


1 


I r) 




/i(0 I r) 


(2.3) 




ViT) I I - ^2 I 't) 


- Z2\t) 


(2.4) 

where we have already used the modular symmetry to simplify the expression. This result 
could be derived by using the Ward identity for the twist operators [T7| . 

For classical part, we need to find all of the classical solutions and calculate their action. 
For convenience we redefine the fields as 

n—1 

X^^’^\z,z) = ^e^^'^X^^\z,z), (2.5) 

i=o 

where 0 < fc < n. For each redefined field X^^'^\z^z), when the argument goes around the 
branch point zi or Z 2 , it gets an extra phase or . Moreover, the boundary condition 


(2.3) changes to 




( 2 . 6 ) 
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where is a vector in the lattice defined by 


n—1 


Ak — 


27ri?^^e " +imj^2), £ Z 

j=0 


(2.7) 


These boundary conditions induce the following monodromy conditions 

i + (f = uW 

J la 

I dzdX^^’’^\z) + (f dzdX^^’^^ = 

•J la J la 

where y^’s are the two cycles of worldsheet torus. In terms of cut differentials 

UJ^^\z) = 'dl{z — Zl I T)~^^~l'dl{z — Z 2 I T)~^'dl{z — (1- )zi - Z 2 I t) 

n n 

Uj 2 ^\z) = di(z - Zl I r)“"'i?i(z - Z 2 I t)~'^^~1'&i{z - -zi - (1 - -)z 2 I r), 

n n 

the classical solutions can be written as 

dxm = dxm ^ 


( 2 . 8 ) 


(2.9) 


( 2 . 10 ) 


Solving the monodromy condition, we get 

(,) _ 

det 

Wl^'^^V2-W^^''^Vi 


a' ' = 


det 


hO^) = -w^^'^\i + wl^^\2 

detIT>) 

det W(^) 


( 2 . 11 ) 


where the W functions are defined to be the integral of the cut differentials along different 
cycles. The definition and properties of the IT functions could be found in the Appendix. 
With these results, the classical action for is just 

1 1 


S{k) — 


Anna' I tT]''^^^lT2 


-(I 1^1 IT^^^^ 1^ + I |2| |2). (2.12) 


The above discussion is the same as the one in |13|. However, in 13 , the classical action 


is further simplified using the relation between W functions. The relation turns out to be 


problematic, as shown in 14 . Furthermore, the lattice translation vi and V 2 were determined 
in 


13 along the way in |9|. But notice that we are considering n-sheeted torus, which is very 


different from the n-sheeted Riemann surface got from two interval complex plane. Actually, 
the translation vectors are much simpler in the n-sheeted torus case. On this Riemann surface, 
the spacial cycles and time cycles on n replicas build the canonical cycles, and all of the cycles 
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on the Riemann surface can be generated by these cycles. Once we fix the monodromy around 
the canonical cycles, we can get all of the monodromy on the Riemann surface. 

For the n-sheeted torus, the monodromy condition can be fixed as 


AiX = 2'KRmj 
A2X = 2Tr Ruj, 


(2.13) 


where 


(1) I • (2) 

rrij = m- + im) 

J J 


? + 


rij = n\ + in 


(2.14) 


are complex integers. If we transform into basis, then 


i=o 

n—1 

= 2ttR " '1 

j=o 


n—1 


e " nij 


Oj. 


(2.15) 


Taking these relations into (2.12) and summing over all the twist fields, we get 

2{k) n—1 n—1 


= E ^(E 


na' 

r=l,2 k 


w. 


W, 


i(k) 

1 


yy cos 27r(j — j )— ■ 


1=0 /=0 


n 


+E 


1(0 




2{k) 


n—1 n—1 J 

E E ~ 


1=0/=0 


n ^ / 


J 


(2.16) 


For convenience, let us define 


n—1 


^n' = E 


k=0 

n—1 


w. 


2(0 


IT 1 


1(0 


B.. = y — 1 

" h w: 


1(0 


2(k) 


I k 

cos27r(j - j )- 
n 


cos 27r(j - j')-. 

n 


We can diagnose the two matrices with the same matrix 


(2.17) 




(2.18) 


and they have different eigenvalues 


U ■ A - U = diag{n 


W. 


2 ( 0 ) 




1 ( 0 ) 


\,n\ 


,wy 


1 




1 ( 1 ) 


1 




1 




yv yy 

U-^ -B-U = diag{n\y^\,n\—y,\,...n\ 

wy wy wy ^ 


(2.19) 
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so 


B-^ = -^A. 

We can calculate the classical contribution of the partition function 


( 2 . 20 ) 


^classical — ^ ^ 


ttR 


n—1 


exp[ 


- w 


(r) (r) 

^ ,ni ^ 

3 ’ 3 


a'n 


r=l,2 fc=0 




2{fc) 




l(fc) 


n—1 n—1 1 

yy cos27r(j — j ) — ■ m'prn!^^) 


i=o/=o 


n ^ i 


ttR 


n—1 




an 


r=l,2 fc=0 




i(fc) 


1 


2(fc) 


vri? 


W 2 
2 


n—1 n—1 j 

^ ^ cos27r(j - j')- ■ nf^n^^i 


i=o/=o 


n 


a n ^ 


k=0 


m 


2(k) 


w, 


l(k) 


n—1 n—1 

EE cos 27r(j 

i=o/=o 


/ k 

- J ) - ■ mjm.i 
n J 


ttR 


rt 

'r). 


an 


k=0 


Wj 

w; 


i(k) 

1 _ 

2(k) 


a 


i?2r 


n 

A;=0 


/ n—1 

an 


R^n 


E 

k=0 


W. 


2(k) 


w, 


l(k) 


n—1 n—1 , 

E E cos27r(j -/)- • njn.>]f 

i=o/=o 

n—1 n—1 

EE cos 2n{j — j 

i=o j'=o 
\ \ 2 

n—1 n—1 J 

y^y]]cos27r(j-/)--pjP^. 

i=o/=o 



1 \ ^ p'^n '' 

nR^^^ 



2 ^ exp < 

\mj,pj 

a'n ^ 
fc=0 



- .7 ) - • m,m ./ 
n 


( 2 . 21 ) 


For the last equation, we have used higher dimensional Poisson re-summation. It is remarkable 
that the classical partition function depends explicitly on the W functions and then on the 
interval. This is very different from the result in |13| . Combined with the quantum part, the 
full partition function reads 


— ^classical ^quantum 


- C 71 


n—1 

n 


ViT) I |2 - Z2\t) 


( ^' 1(0 I ( _ ^' 1(0 I ) ^(1-^) 


I \ ^ p'^n ^ 


^2(fc) 

2 ^ exp < 

\mj,pj 

a'n ^ 
k=0 



n—1 n—1 

EE cos 27r(j 

7=0/=0 


'■i9i(2;i - Z2\t)' 
I k 

- J ) - ■ mjm.i 
n J 


a vr 


n—1 


i?2 


n 


E 


k=0 


2{k) 


w, 


l(k) 


n—1 n—1 

^^cos27r(j -J )- -pjp.i ^ I , 

7=0/=o 


k 

n 


( 2 . 22 ) 


where other coefficients have been absorbed into Cn- This is the main result of this paper. 
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3 Low temperature limit 


For the single interval Renyi entropy at finite temperature, there is an universal thermal 
correction coming from the lowest excitation [Is] 

2A 


6S„ = 


gn 

1 — n 


sin ^ 


resin 


ttI 
nL . 


_ 


(3.1) 


where A is the scaling dimension of the excitations and g is their degeneracy. This relation has 
been checked to be true for the vacuum module in the context of AdS 3 /CFT 2 correspondence 
However, it has shown to break down for the noncompact free scalar, which has a 
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continuous spectrum. In this section, we study the low temperature limit of the partition 
function and Renyi entropy of compact free scalar, and check this relation. 

For simplicity, we assume ^ < 1, that means the momentum mode has lower dimension 
than the descendant of vacuum module. We expand the results with respect to q = 


ry(T) = q^{l + 0{q)) 
i?^(0) = 27rgi(l + 0(g)) 

'di{zi-Z 2 \T) = 2qs sirnr {zi - Z 2 ) {I + 0{q)) 


(3.2) 

(3.3) 

(3.4) 




m 




i(fc) 


= i/3 + 


2i(sin^(t-zi))(sin7r(l- |)(t-zi)) 


dt- 


'21 


sin7r(t — zi) 


+ 0{q) 


(3.5) 


For the summation in (2.22), the leading and next leading contribution come from ruj = pj = 0 
and rrij = 0,pi = ±l,Pj = 0 for j i. When all of rrij and pj are zero, then the exponential 
is 1. When only pi = ±1, the terms on the exponential contribute 


/ n—1 

a 7T 


K^n 


E 

k 

n-1 


Wr 


m 


w, 


i(k) 


f IL—L « 

SE(/^+/ 


- ^ 1 ) smvr(l zi) 


zi sin7r(t-2;i) 

resin ^( 2:2 - zi) 


+ 0(q) 


sin 71(22 — 2 : 1 ) 


• + 0{q). 


Thus the leading and subleading contributions in the summation give 

1 , O sin7r(22 - ^;i) 

l + 2 ree (——577 --)h^ + 0 (e ) 

resin 422 - 21 ) 

and the partition function is approximately 


7 - ^ i ^ 

— Cji n \ . , 

qe sinTTt 


i„(i_ 1),^ . -EslJ. , sinvr/ 

3 ' n.2 ' I 


(1 + 4ree ( 


re sin -I 


) 1 R + 0 (eHi^)) 


(3.6) 


(3.7) 


(3.8) 
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and 


Sn — Cji -|- 


n + 1 


1 


-4n( 


sinvr/ 


3n n — 1 'n sin -I 

n 


/2 

Q 

R e 


'KCip 2'KCX.p 

+ 0{e~ ). 


(3.9) 


This result is in consistent with the universal behavior suggested in |15| . For the complex 
scalar, the central charge is 2. There are four degeneracies for the lowest excitation states, 
with the conformal dimension ( 4 ^, -§^)- 


4 Relation between thermal entropy and entanglement entropy 


In 11 , it was suggested that when the interval becomes very large, the entanglement entropy 


and thermal entropy could be related by 


Sth = lim(S'^;£;(l - e) - See{^))- 

e-5-O 


(4.1) 


In 14,16 , this relation has been proved for a general CFT with discrete spectrum. For the 


free compact scalar case at hand, it has discrete spectrum so should satisfy the relation. We 
can prove it directly by expanding the W functions. 

Since it is not obvious on how to take n —?■ 1 limit on the Renyi entropy, we do not have 
explicit form of the entanglement entropy. Instead, we first try to study the e —)• 1 limit 
and then take the n —)■ 1 limit. We assume that taking the limits of n —)• 1 and e —>■ 1 is 
commutable. 


For small interval, we have 




i(fc) 


= 1 + 0{zi - Z 2 ), 


w. 


2(k) 


= i/3-f- 0(zi - Z 2 ), 


(4.2) 


and thus 


^n. — 






1 


h(T)l^” 

[( 


( exp[-^m 2 - ^ 


y,pj 




^ J^exp[-^m 2 - 


m,p 


a' 


— Cnl 




^>[Z^ + Oil)] 


(4.3) 


For large interval, we first analyze the summation terms in the exponential. We only 
extract the leading terms with respect to zi — Z 2 - Since the eigenvalues for the matrix A 
\,k = 0, ■ ■ ■ n — 1, in the limit zi —)• Z 2 , all the terms are suppressed but the 


are n 


I w. 


2(fe) 


l(fe) 


terms with mi = m2 = ...mn = m,pi = P2 = ■■■Pn = P- This fact is consistent with the 


observation in 14,16 that the excitations on different replicas should be the same in order to 


give nonvanishing contributions in the large interval limit. Therefore the summation yields 

72 ^ R a'n/? 


E 'KR‘^nf3 2 

exp[- —m 

m,p 


(4.4) 
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In this case, when zi ^ Z 2 ^ the prefactor before the exponential goes to 


Cnl 


1 TT, 2sm^r?(r)3 


This factor could be simplified more 

n—l n nk 


1 


j^,2sinf77(i/3)3,2 




n—l , . n—l 

n(2sin-)={/S-i,,(r))2”-«f] 

k=l ^ ^ k=l 


13 he 


n—l 






fe=i 


n ^ r]{inl3)^ 


As 


1 TT I'o • 1 

^||( 2 sm—) = 1 , 

±.L fi 

k=l 


we find that for the large interval 


Zn, — Cnl ® 




(Zi[n/3]+0(C)), 


(4.5) 


(4.6) 


(4.7) 


(4.8) 


where A < 1 and 0{l^) terms coming from the contributions of the primary fields in the 
operator product expansion of two twist operators. 


Now we find the similar results as the ones in 14, 16j. It is easy and straightforward to 


prove the relation (4.1) between the thermal entropy and the entanglement entropy. Actually, 

lim(S'E£;(l — 1) — See{1)) 

1^0 


= — lim 


>1 n — 1 


(logZi[n^] - nlogZi[/3]) 




= Sth. 




Acknowledgments 

The work was in part supported by NSFC Grant No. 11275010, No. 11335012 and No. 11325522. 


10 



















r- 

i ^ i 



\ 



J --^ 

-- 1 .__ 





_e f 1 — 

*2 *1 * 

(a) A cycle 


-«2 • 

(b) B cycle 

1 £2 1 


Figure 1: There are two contour integral paths. 


Appendix A W functions 

In this Appendix, we list some useful results for the W functions and study their properties 
in various limits. The W functions are defined as the line integral of the cut differentials 

— Z 2 I t)~^'i!)i(z — (1 - )zi - Z 2 I t) 

n n 

t)~^^i(z — Z 2 I T)~^^~^'^'dl{z - Zl — (1- )Z 2 I r) 

n n 

— Z 2 I t)~^'&i{z — (1- )zi - Z 2 I t) 

n n 

— Z2 I — —Zl — (1 — —)^2 | '^)-(A.l) 

n n 

As we have chosen the modular parameter to be pure imaginary r = i/3, the W functions are 
related by 

Wl* = = W^, W^* = -W^ = IT|. (A.2) 


= / dz'&iiz-zi 


= / dzdi{z-zi 

Jo 

= [ dzdiiz-zi 
= I dzdi{z-zi 
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Appendix A.l Low temperature expansion 

The low temperature expansion for is 


W, 


i(k) 


= {[ +[ + [ )du'di{z - zi \ t) ^'^'di{z - Z 2 \ t) ^-diiz - {1 --)zi - -Z 2 \ t) 

Jo Jz, Jzo n n 


'0 ^ Zi j Z2 

[ dz'&i(z — zi I T)~^^~"\'di(z — {z 2 — 1 ) I T)e^^)~^'di(z — (1 — —)zi — —Z 2 ) I r) 
Iz 2 -i n n 


_ 7vik 

e 


k. 


k 


2 ^ik Y dz{s\n.'!i{z - zi)) »^(sin 7 r( 2 ; - ( 2:2 - 1 ))) » sin7r(2: - (1 --) 2 ;i -- 2 : 2 ) + 0 (g) 
1 — e » J A ^ ^ 


' /c 

/I dU _k : N-*/ (1--) 

- krzkif " (u-ui) ^ ^\u-U 2 ) ^{u-u\ " u^) 

+ /oo 27rzu 


2TTiu 


—— U 2 )~"{u — u)" + 0{q) 




— Tzi — 
6 ^ 


— 1 + e J(xi 

= l + 0{q). 


J 00 


(A.3) 


In the forth equation we have used a conformal transformation 


u = e 


^2'Kiz 


(A.4) 


and in the fifth equation we have changed the integral contour to the one surrounding the 
infinity and the origin. For we have the expansion 


VF. 


m 


u k \ A* A* 

= / dz‘d{z — zi I T)~"'d{z — Z2 I ~^^'d{z — {—zi + (1 - )z2) \ t) 

Jo n n 

fO f, f, 

= / dz{{l - e 2 ^* 0 -^i))(l - _ q^-2^i{z-z2)^y(i--) 

Jo 

We can analytically expand the equations with respect to and the summations in the 
expansions still converge. There appear several kinds of terms in the expansion 


1 27vmiz m —27rmiz m 27ri{—m-\-n)z 

1, c ^ q t , y 6 


(A. 6 ) 


After being integrated, the first term gives f/3. The second and third terms give 


rO 


^2'Kmiz 


^ „2'Kmiz I i/3_ ^ 


2TTmi 


° 2Trmi 




(A.7) 
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rip 


1 


I ^—2nmiz _ 

Iq —2TTmi 


,m —2iTmiz \iP _ 


1 


-q e 


lo 


—27: mi 


rq^iq-^ - 1 ), 


rip 


^m^2nii-m+n)z ^ 


(A. 8 ) 


(A.9) 


If we are interested in the leading contribution with respect to q, we only need to consider the 
first three kinds of terms. In the end, we have 


W. 


2(k) 


Jo 


+ dz[{l - ^ 

Jo 

/ 200 7 7 

dz[{sm.7:{z — 2 : 1 ))“" (sin 7 r( 2 ; — 2 : 2 ))~^^~"Hsinvr( 2 ;- zi — (1 - )z 2 )) — 1 ] + 0{q) 

-ioo ^ ^ 

(A.IO) 


To deal with the integral in the last relation, we define 

00 7 7 

dz[{smTT{z - zi))~^ (sin 7 r( 2 ; - Z 2 ))~’'^~"Hsinvr( 2 ;- zi - (1 -) 2 ; 2 )) - 1 ], 

-ioo ^ ^ 

(A.ll) 

which only depends on zi — Z 2 - Considering 


dz^F{zi,Z2) 


k 


= TT- 


sin 7 r(l - ^){zi - Z2) n J_i 


(sin7r(2; — 21 )) » ^(sin7r(2: — 22 )) ("^-12) 


d 


dziF{zi,Z2) 


Z2 


■sin7r(l-^)(2i-22) 
we find 


u h h k 

= 7r^ — (1 -) / (sin7r(2;—2:i))“"“^(sin7r(2;— 22 ))”^^””^ cos7r(2;— 22 ), 

^ ^ J—ioo 

(A.13) 


dz^F = - 


1 — V . vr/c 


sin — {z 2 — 2 ^ 1 ) cos 7 r (22 — -^i) — 


dz^F 


n 


n 

—sin —(22 - 2i)sin7r(22 - Zi)dz 2 {- 


sin7r(l-f)(2i-22) 


n ^ 

dz,F 


TT- n 


sin 7 r(l- f)(2i -22)' 


Define 


T = dz^F{zi,Z2), 


(A.14) 

(A.15) 


we have the equation 

dz^T A: cos 7 r ^(22 - 21 ) n - A: cos 7 r(l - ^)(22 - 21 ) cos 7 r( 22 - 2 i) 

—?FI— ~ - h -1" - k -- 7 -2) 

T n sinTT^ n sin 7 r(l - |)(22 - 21 ) sin 7 r( 22 - 2 i) 

which has the solution 

sin^( 22 - 2 i)sin 7 r(l - |)( 22 - 21 ) 
sin7r(22 — 2i) 


(A.16) 


(A.17) 
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Comparing with the direct evaluation of dz 2 F{zi, Z 2 ) for ^2 —)• -Zi, we have 


C = 2i. 


(A.18) 


Therefore, 


and 


F{zi,Z 2) = f 

J Zi 


^^ 2isinf (t - zi) Sin7r(l - ^)(t - zi) 


=i^ + j dt 


sin7r(t — zi) 
^2 2isin^(t-zi)sin7r(l- 


sin7r(t — zi) 

Appendix A.2 Small and large interval limits 


+ 0{q). 


(A.19) 

(A.20) 


To study the relation between thermal entropy and entanglement entropy, we need to study 
the small interval limit and large interval limit of the W functions. For small interval, it is 
simple to calculate. When zi = Z 2 -, the integrand is 1, so 




i(fc) 


= 1 = iP 


(A.21) 


For large interval, we set zi —)• 0, 2:2 —)• 1, 


W, 


i{k) _ 


_ irik 

e ^ 


2 -^ik r dz'&i{z - 2 ;i I r) '^^'di{z - (z 2 - 1) | r) "-i?i( 2 ; - (1 - -)zi - -Z 2 | r) 
1 — e " ■'A ^ ^ 


_ TTZk 

e ^ 


k 


r dzMz) di(z --\r) + 0{zi - (22 - 1 )) 

1 — e n J A ^ 

= I t )+ 0(21 - (22 - !))■ 

2sm7r- n 

n 

Let us evaluate the most singular term in IF 2 when 21 —)■ 0 and 22 —> 1 


(A.22) 


= e [ dz'di(z — zi\t) "'di{z — (z 2 — 1)\t) — —zi — (1 — —) 2 ; 2 |r). 

Jo n n 


,2(k) 


There is a (2 — zi) term in J^i{z — 2 i|r), and a (z — (z 2 — 1)) term in i9(z — (z 2 — 1)|t). The 
most singular terms in the integral come from the integral range near the origin. In this limit 


'di{z — zi\t) ~ 2'it{z — zi)rj{T) 
di{z — {z2 — 1 )\t) ~ 27r(z — (22 — 1))??(t) 
■diiz --zi - {1 --)z 2 \t) ~ --)|r), 


n 


n 


n 


(A.23) 

(A.24) 

(A.25) 


then 


VFo~ e / dz—{z — zi) ri (^z — (z 2 — 1)) (1 — —) | r),(A.26) 

J-iM 27r n 


f^iM 


1 


k. 
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where {z 2 — z\) <€. M (3. Since 

riM 


/ IM 

dz{z - zi)~^{z - {Z 2 - --(1 - (zi - {z 2 - 1)), 

-iM 

2{k) 

the leading singular term in ld^ 2 ' is 


(A.27) 


1^2 ~ i 


sin — k 

-^7?(t)"^T?i(-(1 - -)|t) log {Zl - {Z 2 - 1)), 

vr n 


and 


W. 


2(k) 


w, 


l(k) 


2 sin — 

~ -^(-log(zi - (Z2 - 1))). 

TT 


(A.28) 

(A.29) 
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